Feshbach resonance induced shock waves in Bose-Einstein condensates 
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We propose a method for generating shock waves in Bose-Einstein condensates by rapidly increas- 
ing the value of the nonlinear coefficient using Feshbach resonances. We show that in a cigar-shaped 
condensate there exist primary (transverse) and secondary (longitudinal) shock waves. We ana- 
lyze how the shocks are generated in multidimensional scenarios and describe the phenomenology 
associated with the phenomenon. 
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Shock waves have been thoroughly investigated during 
the last century in several branches of physics including 
fluid mechanics, plasma physics, astrophysics, optics and 
solid state physics (see e.g. Very recently it has 

been recognized that there exists one more type of media 
- ultracold dilute alkali gases in a condensed state, i.e. 
Bose-Einstein condensates (BEC's) - where existence of 
shock waves is possible Isl 0, 0, Q . 

In general, BEC's represent a universal laboratory for 
nonlinear phenomena. In the last few years there have 
been theoretical predictions and experimental realiza- 
tions of many nonlinear structures in BEC's such as vor- 
tices, vortex lattices and vortex rings 0, bright |^ and 
dark yisolitons, collapsing waves 10] and stabilized soli- 
tons to cite a few examples. 

Theoretical studies of shock waves in BEC's so far have 
been restricted to effectively one-dimensional (ID) mod- 
els and explored regimes of free expansion of the conden- 
sate. Moreover, the initial conditions considered in some 
of the previous studies seem to be not experimentally 
feasible. In this Letter we propose and investigate a new 
mechanism for the generation of matter shock waves in 
BEC's which could be easy to implement experimentally. 
As a key ingredient of our analysis we keep the multidi- 
mensional character of the condensate, something which 
IS missmg m previous works |iy,y,lsll2|. 

Estimates for shock wave generation. - In order a shock 
could develop from an initially smooth pulse, the quan- 
tum pressure must be negligible at the initial stages of 
evolution. In this situation the hydrodynamical approach 
(see e.g. 12j) holds and exploiting the analogy one can 
expect the appearance of a breaking point and the sub- 
sequent development of a shock wave in a BEC. In other 
words, for the breaking point to be reached one has to 
require |(V^v^)/v^| <C STrasU where n is the density 
of the condensate and as is the s-wave scattering length 
(considered positive hereafter). The quantum pressure 
can be estimated as 1/Lq where Lq is a characteristic 
scale of the condensate to be estimated as the minimal 



scale between all the characteristic spatial scales of the 
order parameter and its first derivative. This scale must 
be larger than the healing length, ^ = (STrna^)"^/^. The 
hydrodynamic approximation and thus the previous esti- 
mates fail when discontinuities appear since when a shock 
starts to develop Lq ^ because of growth of the first 
derivative near the edge of the shock (although the typi- 
cal size of the condensate is preserved). Due to the role 
of the kinetic energy near the shock edge one should not 
expect formation of real discontinuities (like ones ana- 
lyzed, for example in Q) in a matter- wave system. This 
is similar to what it happens in real fluids where the exis- 
tence of dissipation smoothes out the shock waves. When 
the characteristic size of the wave front becomes of order 
of the healing length, thus making the dispersive term 
relevant, we expect the appearance of oscillations of the 
wave front near the quasi-discontinuity. This behavior, 
which differs from the fluid case, is due to the essential 
differences between the regularization mechanisms of flu- 
ids (diffusion) and matter waves (dispersion). 

The condition Llagn 1 is not sufficient for the gener- 
ation of shock waves since there are two essential factors 
not accounted for by this estimate. First, the density, and 
thus the sound velocity rapidly decrease during the con- 
densate expansion. Second, the largest spatial derivative 
of the condensate is not achieved where the sound veloc- 
ity has its maximum. As an example in the well known 
Thomas- Fermi approximation this point is reached at the 
minimum of the sound velocity. This requires special, 
with finite support 2] or non-smooth 0, initial states of 
the condensate, in order to produce shock waves. 

Our shock- wave generation method is based on the use 
of Feshbach resonances (FR) to control the nonlinear 
interactions. The idea is to make increase with time 
thus increasing the time domain in which the quantum 
pressure is negligible. In order to make qualitative esti- 
mates of the effects expected in this scenario we assume 
that initially the product asfimax is negligible so that the 
initial spatial distribution of the atoms is well approx- 
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FIG. 1: (a,b) Initial (t — 0, thin lines) and final (t — t/, 
thick lines) profiles of a quasi- ID condensate for v — 0.05 
(tf = 0.45) and z/ = 1 (t/ = 0.21). (c) Dependence of the 
breaking time, t^, of the ID shock wave on the aspect ratio v 
for To = 0.16 and tq = 0.02. 
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imated by a gaussian distribution. For simplicity, we 
consider a cylindrically symmetric trap, with transver- 
sal and longitudinal frequencies uo^ and respectively. 
Then the characteristic spatial scales of the problem can 
be identified as the healing length ^, and the transverse 
L±_ = ^/h/muJ±_ and longitudinal Ly = ^yh/mij\\ linear 
oscillator lengths. 

As the s-wave scattering length is ramped up to a larger 
value we get an expansion of the initial state with a ve- 
locity gradient strengthened by the parabolic potential. 
If the characteristic time, tq, of the growth of is small 
enough and the increase of is large enough then at the 
initial stages one can consider free expansion of the con- 
densate. In our case there are two characteristic velocities 
of the expansion were a stands for ± 

or II depending which direction is considered. There are 
four characteristic times in this problem: the two break- 
ing times ta = La/cmax (i-G. the times at which the 
breaking would occur if expansion would be absent) and 
the expansion times = L^/va- We notice that these 
estimates involve Cmax which is the maximal sound veloc- 
ity after the growth of the scattering length is terminated. 
In order the breaking point to be reached experimentally 
in a given direction one should have ta <C Tq, what can 
be rewritten as the condition = L^^asTimax ^ 1 where 
^max is the EEC's density in the center of the trap after 
the increase of the scattering length ^5|. At this point 
it is important to emphasize that due to existence of the 
trapping potential only the regions of large density un- 
dergo a quasi-free expansion while the expansion of the 
low-density regions located near to the periphery of the 
atomic cloud is suppressed by the potential. This is in 
a sharp contrast to the truly free expansion considered 
elsewere |^ H, ^ IE 1^ • Our estimates suggest that there 



FIG. 2: Results of 2D simulations for tq = 0.02, u = 0.05, 
and grid parameters At = 0.001, iVo. = 2048, A^;, = 1024 
on {x,y) e [-30,30] x [-60,60]. (a) Short and (b) long 
time evolution of the derivatives Sx{t) = maxcc,z \dijj/dx\. 
(c) Short and (d) long time evolution of the derivatives 
Sz{t) = maxcc,z \dijj/dz\. Dotted lines correspond to To = 0.16 
while solid lines are for To = 0.14. 



should exist two types of shock waves in a cigar-shaped 
condensate, where <^ ^||- A first shock wave would be 
developed in the transverse direction while a second shock 
wave is expected to arise in the longitudinal direction. 

The model.- Our model is the usual time-dependent 
Gross-Pitaevskii (GP) equation 
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The adimensional quantities used are related with the 
physical spatial (r^) and temporal (r) variables through 
~ ^jl ^ — ^±^- The new wavefunction is related 
to the physical one ^ by V^(x, t) = ^(r, r)^ L\/N ^ where 
N is the number of particles in the condensate. Finally 
U = 4:7rNas/L±^ ag being the scattering length of the 
cold collisions within the condensate and v = uj\\/uj±. 
The normalization for ip is J \ip\'^d^x = 1. 

We ramp up the scattering length according to 



U{t) = 



a 



a [e 



,To/to 



t<To, 
t>To. 



(2) 



In what follows we choose a = 1 without loss of general- 
ity and study the dependence of the formation of shock 
waves on the parameters Tq, tq, i'. 
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FIG. 3: Formation of the first shock along the x direction for 
parameters To = 0.16, ro = 0.02, z/ = 0.05. (a) for 
t = and t = 0.26 just after the first shock has formed, (b) 
Detail of the region x E [1,3] where the irregularities in which 
the shock degenerates are more evident, (c-d) Surface plots 
of \il;{x,z)\ for (c) t = and (d) t = 0.26. In all cases the 
arrows mark the location of the shock wave. 



One- dimensional shocks.- We first consider the cross- 
section of the condensate along one of the directions (ra- 
dial or axial) as an effectively ID case. By solving nu- 
merically the ID version of the GPE (^) we have found 
the dependence of the time of the shock wave formation 
tb on the size of the condensate characterized by the trap 
asymmetry. As it clear from FigC]the shock waves de- 
velop faster in time in systems with smaller aspect ratio 
in agreement with our previous qualitative arguments. 

Two-dimensional shocks.- We have studied the forma- 
tion of shock waves from gaussian ground states under 
variation of the parameters tq, and Tq. Results of typ- 
ical numerical simulations of Eq. are summarized in 
Figs. EE and H corresponding to To = 0.14, tq = 0.02 
thus /7max = U{To) ^ 1100 [Fig. El solid lines] and 
To = 0.16, Tq = 0.02, /7max 3 X 10^ [Fig. Ela,c) dotted 
lines and Figs. IHland^. The increase of the nonlinearity 
leads to an expansion of the condensate which is faster 
along the transverse direction due to the higher compres- 
sion and thus the increased non-stationarity along this 
direction (notice that this is the case t± <C t\\ qualita- 
tively described above, see also 3]). On long time scales 
the effect of the potential is to confine the wavepacket 
thus generating recurrent oscillations. 

As it has been predicted above qualitatively and on the 
basis of the numerical simulations of the ID model, the 
faster transversal expansion implies that the shock wave 
is first generated along the transverse direction and not 
along the longitudinal one as one could naively think be- 
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FIG. 4: Formation of the first shock along the z direction for 
parameters To = 0.16, ro = 0.02, zy = 0.05. (a) |'0(x,^)| for 
t = 0.8 and t = 1.3 after the shock has formed. (c-d) Surface 
plots of \ip{x,z)\ for (c) t = 0.8 and (d) t = 1.3. In all cases 
the arrows mark the location of the shock wave. 




FIG. 5: Surface plots of \^{x,z)\ for (a) t = 2.7, (b) t = 3.0 
and (c) t — 3.5 showing the formation of strong shocks for 
long times due to the oscillations of the wavepacket width in 
the harmonic potential. 

cause of the geometry of the condensate. The formation 
of this shock is seen as an increase of the maximum value 
of the spatial derivatives along x [Fig. Ela)]. This phe- 
nomenon is accompanied by a small sympathetic increase 
on the derivatives along z which does not correspond to a 
true shock. The true longitudinal shock appears later due 
to the slower expansion along this direction [Fig. El^c)]. 

The structures of the longitudinal and transverse 
shocks are shown in Figs. Inland 01 It can be seen how 
the longitudinal shocks lead to a very oscillating behavior 
near the shock edge while for the transverse shocks the 
behavior is smoother, although small amplitude oscilla- 
tions are also present. 

For longer times and due to the effect of the trap- 
ping potential we observe the recurrent formation of very 
strong shocks (see the peaks in Sx-) Sz for t — tt, 27r in 
Fig. El- III Fig- El we observe that the wave becomes 
filamented due to the strongly compression induced by 
the combined effect of nonlinear forces and the harmonic 
potential, a phenomenon which could be used as a very 
efficient matter wave compressor. 

The amplitude of the shocks increases as the nonlin- 
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FIG. 6: Amplitude and breaking time of the longitudinal and 
transverse shocks as a function of [/max for different spatial 
asymmetries. Ax — Sx{tb,x) and Az = Sz{tb,z) with Sx and 
Sz as in Fig. 2. 



earity is increased. In Fig. El we study the dependence 
of the shocks amplitude and breaking times as a function 
of the trap asymmetry u. It can be seen that as the trap 
becomes more symmetric the breaking time for the trans- 
verse and longitudinal shocks become more similar. This 
is a reasonable result since in the case u = 1 the shocks 
should develop at the same time along all the spatial di- 
rections starting from symmetric initial data. We have 
also studied the dependence of the shock formation on the 
rising time of the nonlinearity tq and seen that it plays 
an essential role in the transverse shock formation. For 
instance, taking tq =0.1 and Tq = 0.70, which give the 
same values of /7max as those shown in Figs. I2I3I and 
we have found that there are no transverse shocks (only 
those corresponding to the revival of the oscillation in the 
trap appear). This effect has a simple explanation: since 
Lj_ <C for a definite region of densities Umax one may 
have R\\ ^ R± ~ 1, what means that it is possible to 
reach the breaking point only longitudinally. 

Three-dimensional condensates.- We have also studied 
the formation of shock waves in three-dimensional con- 
densates using the same type of numerical methods on 
grids of up to 256x256x200 points. With this resolu- 
tion we have been able at least to detect the formation 
of the first shock waves. We have found features similar 
to those described for the two-dimensional case with the 
shocks first generating transversely and later longitudi- 
nally if the nonlinearity rises fast enough. 

To conclude, we have proposed a mechanism for gener- 
ating shock waves and described their dynamics. In do- 
ing so we have kept the multidimensional nature of the 
system, something which is missing in previous works 



and could lead to flawed results. We hope that our 
method could guide experimentalists in the field of mat- 
ter waves to observe the formation of shock waves with 
Bose-Einstein condensates. 
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